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Anwendung der Weierstrass’schen Theorie der elliptischen Funktionen 
zur Bestimmung der Bewegung eines materiellen Punktes auf einem Kreise, einer 
Kettenlinie und einer Parabel. 


Vorbemerkungen, 


Während die neuere mathematische Litteratur keinen Mangel an trefflichen Werken über die 
Theorie der elliptischen Funktionen aufweist, fehlte es bisher an einer geeigneten Zusammenstellung 
von Aufgaben zur Anwendung dieser Theorie. Die meisten Anwendungen giebt Schellbach’s grund- 
legendes Werk „Die Lehre von den elliptischen Integ alen und den Thetafunktionen“ (Berlin 1864); 
Durée behandelt in seiner „Theorie der elliptischen Funktionen‘ (4. Auflage Leipzig 1887) das ebene 
und das sphärische Pendel; Thomae untersucht in seiner „Sammlung von Formeln, welche bei An- 
wendung der elliptischen und Rosenhain’schen Funktionen gebraucht werden“ (Halle 1876) die Be- 
wegung eines materiellen Punktes auf einem Kreise und einer Parabel; Hermite unterwirft in dem 
Aufsatze „Sur quelques applications des fonctions elliptiques“ (Comptes Rendus 1877) das Problem 
der Rotation eines Körpers um einen festen Punkt einer erneuten und ausführlichen Betrachtung; 
endlich giebt Bobek in seiner „Einleitung in die Theorie der elliptischen Funktionen“ (Leipzig 1884) 
auch eine Anwendung der elliptischen Funktionen auf die Geometrie, nämlich auf die ebenen Kurven 
dritter Ordnung. Ausserdem liegen noch in den verschiedenen mathematischen Zeitschriften, Programm- 
abhandlungen und Dissertationen Arbeiten vor, welche sich mit der Lösung solcher Aufgaben be- 
schäftieen, die auf elliptische Funktionen führen. 

Erst vor Kurzem ist diesem schon oft empfundenen Mangel an einer einheitlichen Behandlung 
der mittels elliptischer Funktionen lösbaren Probleme durch das Erscheinen des Halphen’schen Werkes 
„Theorie des fonctions elliptiques et de leurs applications“ (Paris, Gauthier-Villars) abgeholfen worden. 
Der erste im Jahre 1886 erschienene Band dieses grossartig angelegten Werkes (Gross-Oktav, VIII. 
und 492 Seiten) enthält die Theorie der elliptischen Funktionen; der Ende vorigen Jahres heraus- 
gegebene zweite Band (#59 Seiten) behandelt die Anwendungen der elliptischen Funktionen auf Probleme 
der Mechanik, der Physik, der Geodäsie, der Geometrie und der Integralreehnung; ein dritter Band 
endlieh soll sich mit der Theorie der Transformation und mit den Anwendungen auf die Algebra und 
Arithmetik beschäftigen, sowie einen historischen Überblick über die elliptischen Funktionen geben. 
Während in den meisten der oben angeführten Abhandlungen und Werke die von Jacobi, dem Be- 
eründer der Theorie, eingeschlagenen Methoden zu Grunde gelegt sind, fusst das Halphen’sche Werk 
durchaus auf den neueren Weierstrass’schen Bezeichnungen und Methoden. 

Es ist dies auch ganz natürlich. Denn, obwohl für manche Probleme, wie z. B. das Problem 
des gewöhnlichen Kreispendels, die Jacobi’sche Methode, für andere Probleme, wie z. B. das der 
Rotation eines schweren Körpers um einen festen Punkt die Weierstrass’sche Methode sich wie von 
selbst darzubieten scheint, so unterliegt es doch keinem Zweifel, dass diese letztere Methode vor jener 
mannigfache und gewichtige Vorzüge besitzt. Ihr wesentlichster Vorzug besteht, wie ich schon in 
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meiner Abhandlung ‚Neue Lösung eines Rotationsproblems‘ (Zeitschrift für Mathematik und Physik, 
Bd. XXVI., pg. 104 ff.) hervorgehoben habe, darin, dass sie die so unerquickliche Rechnung mit 
elliptischen Integralen dritter Gattung vermeidet und überhaupt alle jene Methoden überflüssig 
macht, welche zur Reduktion elliptischer Integrale in reeller Form dienen, sowie die Untersuchung 
der mannigfachen Fälle, die bei dieser Reduktion zu unterscheiden sind. Andererseits ist man bei 
Anwendung der Weierstrass’schen Methode nicht in Verlegenheit um die zweekmässigste Substitution, 
welche das betreffende elliptische Differential auf die sog. Normalform überführt; der einzuschlagende 
Weg ist hier bestimmt vorgezeichnet, und auch hierin liegt ein gewiss nicht zu unterschätzender Vorzug. 

Leider kommt dieser letzte Vorzug in der Halphen’schen Darstellung nicht zur vollen Geltung; 
es fehlt derselben das einheitliche Prinzip, welches den von Herrn Weierstrass in seinen Vorlesungen 
gegebenen Lösungen zu Grunde lag und denselben eine so hervorragende Einfachheit und Klarheit 
verlieh. Die Herausgabe der Weierstrass’schen Vorlesungen über die Anwendungen der elliptischen 
Funktionen würde daher auch heute noch einem dringenden Bedürfnisse abhelfen, umsomehr, als es 
nur bedauert werden kann, dass es einem französischen Mathematiker zuerst vorbehalten war, die von 
unserm grössten deutschen Mathematiker gewonnenen Resultate in allerdings durchaus selbständiger 
Weise zu verarbeiten und der Öffentlichkeit zu übergeben. 

Die vorliegende Abhandlung verfolgt einen dreifachen Zweck: Erstens soll sie an einigen der 
einfachsten mechanischen Probleme, die bereits von anderen Mathematikern wenn auch auf anderem 
Wege gelöst worden sind*), die Weierstrass’sche Theorie zur Anwendung bringen; sie soll ferner nach 
dem Yorgange von Schellbach diese Theorie durch einige vollständig durchgeführte numerische Be- 
rechnungen erläutern, da durch diese das Verständnis der Theorie wesentlich gefördert wird; endlich 
soll sie zeigen, wie bei diesen numerischen Berechnungen die hyperbolischen Funktionen in 
gleicher Weise zur Anwendung kommen, wie die gewöhnlichen Kreisfunktionen. In dieser letzten 
Beziehung wird, wie ich hoffe, auch dem schon so vielfach behandelten Probleme des Kreispendels 
eine neue Seite abgewonnen werden, welche der Aufmerksamkeit der Mathematiker bisher entgangen 
zu sein scheint. Zwar hat bereits Hoüel in einem kurzen aber inhaltsreichen Aufsatze „Notices sur 
les fonctions hyperboliques et sur quelques tables de ces fonctions“ (Nouvelles annales de mathématiques, 
Jahrgang 1864, pg. 416 ff.) auf die Verwendbarkeit der hyperbolischen Funktionen zur Berechnung 
elliptischer Funktionen hingewiesen und an einem numerischen Beispiele zu Jacobi’s bekannter Ab- 
handlung „Sur la rotation d'un corps“ (Jacobi’s gesammelte Werke, herausgegeben von Weierstrass, 
2. Bd. pg. 329) erläutert; ferner hat Günther in seinem Werke „Die Lehre von den gewöhnlichen und 
verallgemeinerten Hyperbelfunktionen‘‘ (Halle a./S. 1881) andere interessante Anwendungen der hyper- 
bolischen Funktionen auf die Theorie der elliptischen Funktionen gegeben. Doch jene ausgedehnte 
Verwendbarkeit, welche die hyperbolischen Funktionen bei den meisten auf elliptische Funktionen 
führenden Problemen finden können, ist auch in diesem an litterarischen Nachweisen überaus reich- 
haltigen Werke nicht hervorgehoben worden, und in Halphen’s oben besprochenem Werke ist von 
hyperbolischen Funktionen überhaupt an keiner Stelle die Rede. 


*) Die Lösung des Pendelproblems befindet sich u. A. in den oben angeführten Werken von Schellbach, Durege und 
Thomae, sowie in dem vorzüglichen Lehrbuche der analytischen Mechanik von Rausenberger (Leipzig 1888); die Bewegung 
eines materiellen Punktes auf einer Parabel behandelt Thomae a. a. O.; endlich hat im Programm der Realschule I. Ordnung 
zu Stralsund, Ostern 1873, Herr Dr. Gentzien die „Bewegung eines materiellen Punktes auf einer gemeinen Kettenlinie 
bestimmt, 


I. Absehnitt. 


Die Bewegung eines materiellen Punktes auf einer Kreislinie. 


Zur Bestimmung der Bewegung eines materiellen Punktes in einer Ebene unter dem alleinigen 
Einflusse der Schwerkraft braucht man nicht auf die Differentialgleichungen der Bewegung zurückzu- 
gehen; es genügt vielmehr das Prinzip der Erhaltung der lebendigen Kraft in Verbindung mit der 
Gleichung der Bahnkurve und den gegebenen Anfangsbedingungen. 

Die Gleichung des Kreises sei unter Zugrundelegung eines rechtwinkligen Koordinatensystems x, z: 

a: Sr 
sodass | den Radius des Kreises oder die Pendellänge bezeichnet und der Anfangspunkt der Koordinaten 
mit dem Mittelpunkte des Kreises zusammenfällt; die z-Axe besitze die Richtung der Schwerkraft, 
die Beschleunigung werde wie üblich mit g, und die Geschwindigkeit des beweglichen Punktes zur 
Zeit t mit v bezeichnet. Dann liefert das oben genannte Prinzip folgende Gleichung: 


'dx\? CEV 
o ve (‘ ) H (F) 三 2gz +h. 
CS Ke l E: Séch 


Zur Bestimmung der Konstanten h werde angenommen, dass der bewegliche Punkt für z = |, 
also beim unteren Durchgange durch die Vertikale, die Geschwindigkeit v, besitze. Die Gleichung (2) 


geht alsdann über in e 
vw = 2gl + h, sodass 


. 8D b = yi agl 
ist. 
Durch Differentiation der Gleichung (1) nach t ergiebt sich 
dx; dz 
X +z 
dt dt 


taj En eye Sin GI 
dt. x dt 1? — z? dt) ; 


substituiert man diesen Ausdruck in Gleichung (2), so erhält man zur Bestimmung von z folgende 
Differentialgleichung I. Ordnung: 


a (Z= (1-3) Ce +h = Ro, 
at E 


in welcher unter R (z) das Polygon 2 (oder An) Grades auf der rechten Seite der Differential- 
gleichung zu verstehen ist, also ein Polynom, welches in allgemeiner Gestalt die Form 

5) R(z) = Azt + 4Bz? + 602? + 4B'z + A’ 
besitzt. Im vorliegenden Falle ist 
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\ g z h ' g ' 
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Die Wurzeln der Gleichung R (z) = 0 sind hier sämtlich reell; ihre Werte sind (in irgend welcher 
Grössenfolge): 


V 
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W= ahh 

5- = 21l ist. Um demnach 
2g < 


die 3 Wurzeln, die wir jetzt mit e, @, gx bezeichnen wollen, ihrer Grösse nach zu ordnen, indem wir 


Die grösste dieser 3 Wurzeln ist o, ferner ist ge, je nachdem 


1) 4 > a > 


annehmen, müssen wir 2 Fälle unterscheiden. Im ersten sei 


y 2 
一 <2] 
Zo d 
dann ist 
vor 
Gel, Gs el, Gezei 
2g 
im zweiten Falle sei 
e 
Vo 21 
CR Sick 
dann wird 
g Vo 
Mash Ban) ray = ， 
28 


Im ersten Falle ist, wie aus den Gesetzen des freien Falls und aus dem Prinzip der Erhaltung 
der lebendigen Kraft hervorgeht, die Bewegung eine oseillatorische; im zweiten Falle durchläuft der 
materielle Punkt die ganze Kreisperipherie und führt eine fortdauernde Rotation aus. Im Grenzfalle 
v? 
2g 


— 2] nähert sich der bewegliche Punkt dem höchsten Punkte der Kreisperipherie asymptotisch. 


Die nächste Aufgabe besteht darin, die Diflerentialgleichung 


EA 本 
CG SE (z) 
时 二 

SE eg ; 


zurückzuführen, in welcher S das Polynom 3%" Grades 


S = 48° — g8 — g, = 4 (s — e) (8 — ea) (S — e) 


auf die Weierstrass’sche Normalform 


bezeichnet. Es geschieht dies, wenn man unter z, eine willkürliche Konstante versteht, mittels der 
Substitution: 
VR (z) VR(z) + Biz , #R’(z) 


8) í z TAR (z) 
SS 2(z 一 2Zo) 2 一 20 = (zo) 
- Go : r dz ds S 
welche das allgemeine Integral der Differentialgleichung 了 = t Ze darstellt *). 
YR (z2) yS 
Wählt man insbesondere für z, eine Wurzel der Gleichung R(z) = 0, so geht die vorige 
Gleichung über in: 
Riis 要 
9) Ee + ( 0) ab 2 R’'(z,) **) 


Z 一 Z0 


*) Die Herleitung obiger Formel findet man in der inhaltsreichen Abhandlung von G. Mittag- Leffler „ièn metod att 
komma i analystik besittning af de elliptiska funktionerna“ (Helsingfors 1876), sowie in der Abhandlung meines hochverebrten 
Lehrers, des Herrn Prof. Scheibner „Zur Reduction elliptischer Integrale in reeller Form“ (Separatabdruck aus dem XII. Bande 
der Abhandlungen der mathematisch-physikalischen Klasse der Kgl. Sächsischen Gesellschaft der Wissenschaften) Artikel 5. 
Vgl. auch Schlömilch’s Compendium der höheren Analysis, 1, Bd. § 109; Halphen Le 2. Bd., 9. Kapitel. 

**) Vgl. Scheibner 1. e, Art. 10. 


q R (zo) 


oder umgekehrt: F 
10) zz Su + DEER CCS 


und dies ist die wichtige, bei jedem mechanischen oder geometrischen Probleme wiederkehrende 
Substitution. 
Wir wählen für z, die Wurzel l, also denjenigen Wurzelwert, den die Variable z sowohl im 
ersten wie auch im zweiten Falle im Laufe der Bewegung anzunehmen imstande ist. Dann ist 
FRO iR () 
ee s— IR”) ` 


Zur Ermittelung der Werte von o, e, €, hat man diejenigen Werte von s zu bestimmen, für 


11) 上 去 了 RD; z=1+ 


welche S verschwindet. Nun folgt aber aus (11) und (4) 


ds IR'() dz | 
= — I als 
dt Be en 
Se (2 = [# R’(1)]? l Cal Ze [ 1 ul. (z +1) (2 gz +h) j 
dt (2—1)* \dt, 41 (z — 1)? 
Dieser Ausdruck verschwindet dann und nur dann, wenn z einen der Werte 
h Vae 
on | eS — 1 SC KC 
eg: 28 
annimmt. Die entsprechenden Werte von s, die einstweilen mit e, e”, e" bezeichnet werden mögen, sind 
1 1 1 > l > or 
Ua er H E, mr KE R ] Na n ER JA r 
e ai H Ui. e za (1) gi (D e SL (1) RE a), 
oder da 
Zi: 
| RO) = — y 
12) a 
Av Ho 
uh 一 u 0 e= > 
=. H l 
ist: d 2 E 9 
一 Vo“ Së 8 A vo E H ida T Vo SS g , 
ra Ee RPT 3I 


Von diesen drei Grössen, die wieder sämtlich reell sind, ist offenbar, da R’ (1) wesentlich negativ 


ist, e die kleinste; da nun p 
vw — 4gl 


41? 


H HI 
e — e 一 -一 
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: : Y ES 2 s i 

ist, so wird e” = e", je nachdem SC = 21 ist, d. h. je nachdem der zweite oder erste Fall vorliegt. 
d £ 
g 


Setzt man daher, wie üblich, die Reihenfolge 
13) o > ea > Ge 
fest, so ist im ersten Fallo ẹ = e”,  =e, Ga =e 
und im zweiten Falle e = e”, e = e”, @& = e’, 
d. h. es ist, wenn wir wie bisher die auf den ersten bezw. zweiten Fall bezüglichen Formeln mit den 


Indices a bezw. b bezeichnen: 


e 2 E ep 7 2 e 2 
ve | 2g y g Yo g 1. 
e = — geng +7, fe zë b = nn E SE UN 
1217 31 z bo ai 121° 31 24 
只 
CS CN H g Vo” 
e — Dn = 一 一 fy se Ba = mr tee ee EE e 
2 3 41?’ 1 d l l - l 4 1? ? 


E ee ER e EN en vm E Er Een 
13») |a 6? 31° Tr pae Ke = Bi S n~ m 
3, 5 


Durch Substitution (11) ist nunmehr die Differentialgleichung (4) übergegangen in 


Ian? 
14) (1) = 46 — 8) — e) (5: — 8), 


deren vollständiges Integral 

15) s = p(t —t,) 
ist, wo p die bekannte Weierstrass’sche doppelt periodische Funktion mit den primitiven Perioden 
2w, 2w und t, eine willkürliche Konstante bezeichnet, welche von den Anfangsbedingungen der 
Bewegung abhängt. Aus den Gleichungen (11), (12) und (15) ergiebt sich somit als das vollständige 
Integral der ursprünglichen Differentialgleichung (4): 


vo? 
有 ai 


p(t to) — e3 
Bei den weiteren Rechnungen werde ich mich auf die „Formeln und Lehrsätze zum Gebrauche 
der elliptischen Funktionen, nach Vorlesungen und Aufzeichnungen des Herrn K. Weierstrass bearbeitet 
und herausgegeben von H. A. Schwarz‘ (Göttingen 1885) beziehen. Die in denselben angewandten 
Thetafunktionen sind definiert durch die Reihen 


L . P . e B - 一 
Ate, ol = 24? (sin væ — q? sin 3 væ + qê sin 5 va F...) 
1 3 7 i 
= (v, q) = 2 q4 (cos væ + q° cos 3 vr + qê cos 5 væ +... 
17) | ma | | -q 
(v, q) = 1 + 2q cos 2 væ + Sai cos 4 væ + 2q” cos 6 va +... 
Ae, ai = 1 — 2 q 0082 vr + 2 q4 cos 4 væ — 2 q? cos 6 vm +... , 
ferner ist 
wni 
i ge 一 日 Ce ei 一 es 
) ee, ee 2 a?ac = 
18 el — Ge e, — Ge 
| e + e& + e = 0, ra e 7 


Zunächst ergiebt sich, da nach der Formelsammlung Art. (18) 


La vm %@) 
on) es = a° (u) 
ist: 
or, a 
21 0° (t— to)’ 
also 
Se ot — to) 
| 21 o (t — to) 


Nun ist i der Cosinus desjenigen Winkels, den das Pendel zur Zeit t mit der z- Axe, d. h. mit 


der Richtung der Schwere bildet; bezeichnet man diesen Winkel mit y, so erhält man 


E a N y ve Vo 2 lt — to) r 
1 cos Y - Zem a gn PETI = to) ; 
mithin 
sin 7 = +3, Al; a R 
2 21 a(t —t) 


Geht man endlich von den g- zu den $-Funktionen über (Formelsammlung Art. 34 und 35) und setzt 


19) u= Z > ‘o 
"so wird Se 
mr — m, J(u) E J(u) _ 1 Vom. 1 Ain) EE 1 9 (u) f 
2 — 70), S, — læ 3(0)9,(0) (u) SH 2 ~ (u) 


oder endlich, da nach (13a) und (13b) in beiden Fällen 


vim 
(Te — es) = 47 
ist: d ; 
20) sin m: Hp Voed J(u) = 
2 dal J(u) 


1. Fall E <2l). 


Aus den Gleichungen (13a) ergiebt sich 
í e; vg” 
2 2 3 0 
a E , also 
ci es 4 gl 
? a “W j J(u) 
20a) sin = e Ta ez 
2 (u) 
Nimmt man an, dass für v = 0 z = z sei und bezeichnet den Winkel, den das Pendel in der 
höchsten Lage mit der positiven z-Axe bildet, (die sog. Amplitude) mit æ, so ist nach Gleichung (2) u. (3) 
h = — 2 gz 
Vo = 2 g l Zo) 
R vo? l — z [ ( =) H i 
pa e = 1 一 = (1 — cosa) = sin’<, 
4 gl 2] 2 ] 2 2 
also 
C 
9]. 二 2 
21a) x sin, - 


Schliesslich ist noch der Wert der Konstanten tọ und das Vorzeichen der Quadratwurzel in (20a) 
zu bestimmen. Rechnet man die Zeit von einem derjenigen Momente an, in denen das Pendel durch 
die Vertikale geht, sodass für t = O auch ui — 0 wird, so ist 

22a) to = 0 
zu setzen. Nimmt man ferner an, dase zu Anfang der Bewegung % positiv sei, so ist in (20a) das 
positive Vorzeichen zu wählen, sodass definitiv 
Q f 
23a) sin S SEN (u) rE 


; 2 Au)’ 2 w 
wird. V 


Da, wie aus den Gleichungen (17) hervorgeht, 


is +1) = — I) 
AY +1) = (y) 
; A $ J(u S 2 < a - a o 
ist, so besitzt der Quotient en als Funktion von u die reelle Periode 2, als Funktion von t hingegen 
ol 


die reelle Periode 4w, d. h. das Pendel befindet sich immer nach Verlauf von 4w Zeiteinheiten wieder 


> 
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in demselben Bewegungszustande, oder m. a. W.: Die ganze Schwingungsdauer ist gleich 4 w. Be- 
zeichnet man die einfache Schwingungsdauer, d. h. die Zeit eines Hin- oder Herganges mit T, 
so ist demnach 

24a) DN 2 w, 


und œw selber stellt die Zeit dar, die das Pendel braucht, um aus der höchsten Lage in die tiefste 
(vertikale) Lage zu gelangen, oder umgekehrt. 

Will man als Anfangspunkt der Zeit nicht wie oben einen Durchgang des Pendels durch die 
Vertikale wählen, sondern einen derjenigen Zeitmomente, in denen das Pendel die höchste Lage besitzt, 
so hat man, wie soeben gezeigt wurde, 

25a) bo w 
zu wählen; man erhält alsdann 


ll 一 (u) t— o t 
sin 二 = x: --—, u= = — 5, 
2 lu) 2 w 20 2 


sodass in der That für t = œw wy = 0 wird. Nun ergiebt sich aber aus den Definitionsgleichungen 
der $-Funktionen, dass 


F (v — è) = — Iy) 
Hv — $) = Hl) 
ist; somit ist bei dieser Konstantenbestimmung 
v g2 t 


26a) sin N 


—a ZS < 5 
2 | J(u)’ 2 w 
Zur Bestimmung der reellen Periode 2 w oder der Schwingungsdauer dient eine der drei Formeln 


(Formelsammlung Art. 55): 


20 /一 E à , 4 as 
一 Ve 一 es v (0) = |2 gi (I1+qQ9?+g°’+gqg"® +...) 
2 w j àa e 9 nå 9 oa 2 

27) = Vë — 8 I (0) = (1 -+ 2q +20 +20 +...) 
2 w 


— Ve ën = Ho) = (1 — 2q + 2q — 20? + ...), 


in denen das Wurzelzeichen positiv zu nehmen ist. Wählt man, um zu dem bekannten Ausdruck für 
die Schwingungsdauer zu gelangen, die mittlere dieser Formeln, so wird mit Rücksicht auf (13a) 
J a 


25) T= 2o S | El = / > -1+2q9q+2g + ...). 


Lë 


Le) 


Für sehr kleine Werte der Amplitude c nähert sich, wie im Verlaufe der numerischen Rechnung 
gezeigt werden wird, q der Grenze 0; man hat also für solche Werte von o die aus der elementaren 


fi 
Mechanik bekannte Formel T = a | . — 


or 
> 


Vorbereitungen zur numerischen Rechnung. 
(Die lineare Transformation der $-Funktionen.) 


Ist der Ausschlagswinkel oder die Amplitude œ gegeben, so ergiebt sich nach Gleichung (21a) 


a: a q 
x = ns: X = ll > cosg und vi —= el. si, 
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Somit sind zunächst die Grössen e, €g, Ge bekannt. Um ferner die reelle Periode 2 w (d. i. 
die einfache Schwingungsdauer) zu berechnen, bedarf es zufolge der Gleichungen (27) erst der Berechnung 
von q. Aus (18) und (27) ergeben sich aber folgende Beziehungen zwischen q einerseits und x, x’ 
andererseits: 


es — êg 4 ( 1 Ek q? / ) 
-这 一 = qË. 一 一 一 -一 
Bon enger 
29) 
— 99 I.9 ci 2 
- /el — es ( 2q + ZO + ==) 
x = I " = ; 
a ıI+2gq+2q4 -+.... 


durch Umkehrung einer dieser beiden Gleichungen, etwa der zweiten, erhält man, wenn zur Abkürzung 


30) 2 Gef e 
= 1+ Val 
gesetzt wird: 
AU ga = Ak 22 + 152° + 15028 Ek OTA +. 


Ist somit q gefunden, so liefert jede der drei Gleichungen (27) unmittelbar den Wert von 2 œw, 
d. h. die Grösse der einfachen Schwingungsdauer des Pendels. 

Der weitere, Verlauf der numerischen Rechnung lehnt sich, wie weiter unten gezeigt werden wird, 
unmittelbar an die Gleichungen (17) und (26a) an, bietet also zunächst keinen Anlass zu weiteren 
allgemeinen Bemerkungen. 

Letzteres ist jedoch der Fall mit der Gleichung (31), die zur Berechnung von q dient. Diese 
Reihenentwickelung konvergiert im allgemeinen sehr gut, besser als sämtliche in der elementaren Ma- 
thematik vorkommende Reihen. Benützt man bei der Rechnung siebenstellige Logarithmentafeln, so 


A 7 Q b . R e 5 ` $ 
kommt für solche Werte von — die kleiner als 38° 15’ sind, schon das 2'° Glied der Reihenent- 


一 


wickelung nicht mehr in Betracht, während, wenn — grösser als der soeben angegebene Winkel aber 
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kleiner als 66° 47’ ist, die beiden ersten Glieder obiger Reihe genügen **). Also nur für verhältnis- 
mässig sehr grosse Werte von « braucht man zur Berechnung von q mehr als die beiden ersten Glieder 
jener Reihe. Es ist nun die Bemerkung von Wichtigkeit, dass man in keinem Falle mehr als die 
beiden ersten Glieder braucht. 

Um diese Behauptung zu beweisen, müssen wir auf eine Umformung der Definitionsgleichungen (17) 
für die $-Funktionen eingehen, die unter dem Namen der „linearen Transformation“ bekannt ist. 

f; 》 
Führt man in den Gleichungen (17) statt der trigonometrischen Funktionen die Exponential- 
5 5 

funktion ein, indem man 


und ausserdem q = e setzt, so gehen dieselben über in 


*) Die einfachste Herleitung dieser Reihe (von den Konvergenzbedingungen abgesehen) giebt Schlömilch in seinem 
Compendium der höheren Analysis, 2. Bd. pag. 438; vgl. auch Schellbach 1. c. § 43, Bobek Le pag. 211 ff., Thomae 
„Elementare Theorie der analytischen Funktionen“ (Halle 1880) pag. 125. 

**) Die oben angegebenen Werte ergeben sich durch eine einfache Rechnung aus den Gleichungen 225 = 0,00000005 
bezw. 1549 一 0,00000005. 


9* 
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l 1 + n (n+ 4)? (2n + I)vai 1 ER n (UnL #20’ + (2n + 1)vo) j! 
AA Or DTA (—1)q e zs e Aal ll e w 
. Op ds Do 
EP? npp (2n + Dest SZ (ntp + (n Leef 
Zeie, dq) = n4 e == RS. zi 
b 一 名 — 0 
32) + CR + 0 mi 
n? 2nvai (n2w -+ 2 n vw) 
zs(V， q) SS nq E - KS © 和 
=Q A 
= n n? 2nvai < n (nzw 十 2nvow) zi 
am = Dada e = J,e o, 
一 Oo — © 


Die weitere Umformung möge auf eine dieser #-Funktionen, etwa auf 3, beschränkt werden. 
Schreibt man die dritte Gleichung in der Form 


vim 4 DC vo j w'ni 
ge ` el 7 $ 
DH e Ee EE w w 
Js (V, 9) > .e H 
on 


und führt eine neue Variable t ein mittels der Gleichung 


vo t w 
33) 一 一 oden y = zdis 
7 Ei ER 
w ZT á DIT 
so ergiebt sich e Ki 
wt? +o wni t wä 
7 $ e 十 d 
s wt t wni wW | SE? j 
$, ok sp S e 2 
"Vë wn gea 
SEA 
oder ; R i 
w t* - +o wni a 
Se Ss { wri 9 ml S , o ( RE: Ss d 
P(t) = e ea oae > a6 - 
207 = 


Diese Funktion gx bleibt ihrem Werte nach unverändert, wenn man n durch n + 1, oder, 
was auf dasselbe hinauskommt, t durch t 十 2 sr ersetzt, besitzt also die Periode 2 æ. Daher kann 
dieselbe in eine Fourier’sche Reihe entwickelt werden von der Form 


CC 
mti 
P(t) m am © S 
al &, 
wo 
Lä 2 — wni t ə 
1 RE — mti l ER ach, E (tal — mti 
t á 
am = ES HOP dt / € dt 
一 区 O — zn 
ist. Um den Summenindex n unter dem Integralzeichen wegzuschaffen, braucht man nur 
t 
SER u 
zu setzen, dann wird Op 
e LL wmi A 
iea en a u? — 2 muzi 
RE J ei du. 
gesellt 
O n— $ 


Die Glieder dieser Summe schliessen sich mit ihren Grenzen stetig an einander, sodass sich am 


in der Form darstellt 
+ u? — 2 muzi 
) 
Au = j e"? du. 


Nun ist 


mithin wird 


und somit 


nm d 


wt? 


os (t) = ei wi s A. ( 


„+ DC 


00 


oder endlich mit Rücksicht auf die Substitutionsgleichung (33) 


Die Summe auf der rechten Seite dieser Gleichung 


wni ,„ 
. y2 


F +% 
| wi 


w 


wri 
r w 
q e 
2 wi 
eg AC 
w) 


setzt, sodass die symmetrischen Relationen bestehen 


so ergiebt sich 


Lé 
> w 
36) wv - 
] 


wri 
D v<“ 
w 
e 1 I (y 
37) e 
wm 
r va 


2 nvy'n 
e 


E Ve af pa mer d a P via 


und Logo, Log o = 


E ee 


/一 (ia)， 


zi 
(— n?w + 2 nyw) 
w 
S d e 
7 n S 


hat wieder die 
man nähmlich neben den Argumenten v und q zwei neue Argumente v’ 


aber, wie aus der Theorie der D Funktionen bekannt ist, 


f b? 
"eg ax? -+ bx Ek jr | e4 a 
d N ` 
a 
j wati 
Hit m? 
E DAR ON 
E: 
w 
See „wzi - 
w't (mi > 二 erg 十 nti 
Yi - / a e 
2 wn l vm u I 


D 


e *), 


und analog 


Form einer #-Reihe. 


Führt 


Be 
und o ein, indem man 


t- 7) fin 
a a E 


*) Es sind hier und im folgenden die Bezeichnungen von Gronau angewendet worden, deren sich dieser um die Theorie 


der hyperbolischen Funktionen hochverdiente Danziger Mathematiker in seinen Schriften über hyperbolische Funktionen bedient, 


insbes. in seinen „Tafeln für sämtliche trigonometrische Funktionen der cyklischen und hyperbolischen Sektoren“ (Danzig 1863, 


im Selbstverlage der Danziger Naturforschenden Gesellschaft). Danach werden die Briggs’schen Logarithmen mit dem Zeichen log, 


die natürlichen mit dem Zeichen Log; die eyklischen Funktionen mit sin, cos ete., die hyperbolischen mit Sin, Cos ete. bezeichnet. 


14 
十 op 
e A n (n 2 (2n + 1)v ia — v'm s dv — vn 
ihai = 2a A-1) q TE Te qile E e )— acle SE ES 
SE? 
F 
7。 Ke Im LE (Gnt (Deia vn — vn CC Ivan 一 8v'r 
RT a E= 2,4 e = qt (e +e ) 十 q (e +e ) +e, 
00 
d. bh, 
gs(viq) 王 1 十 2q' Co82vTr 十 29'0o84Y 六 十 29q20os6v 关 十 .. 
38 ii, ol = 1— 2 qg'Cos 2 v'm + 2 q'*Cos4v'nr— 2q °Cos ô vat... 
38) i9 (v'i, q’) = 2 g'ł(Sin v'æ— q °Sin 3 v'm + q'ĉSin 5 v'a F...) 
RR MN t (Cos v’ -+ q'°Cos 3 v'm + q'Cos5 v'm +...) 
ist. 一 Die Formeln (37) gelten für beliebige reelle oder komplexe Werte von œw und oi: der Wurzel- 
- / è 
Jwi . ; : S SC : è 
grösse | 一 ist, wie bei der Bestimmung des -Integrals in (34) gezeigt wird, derjenige ihrer beiden 
wW j 


£ 
- e e o - w 208 ; S 
Werte beizulegen, dessen reeller Bestandteil positiv ist. Hier, wo w und 一 positiv reell sind, ist 
i ; 


also jene Wurzel ebenfalls positiv zu nehmen. 
Für v = o wird auch vd = o. und demnach 


wW 
wi 
39) be Lët d) = | SC (0,4) 
E: 
WOU 
(0,9) = | a (0,9), 
während 2, (o) identisch verschwindet. Die Gleichungen (27) gehen daher über in: 
2 w we 9 d ) wi 9 2( n 
— Ves — e = 2 (0,0) = —% (0, 0 
5 Fez 3 2 1 ek EC Er S 
2 w = wi 
Ye, 920,0) — ar 2s (0,q ) 
2 w — S DI, 
Ve, — e = 9%0,q) 三 一 92(o,d)， 
w 


aus denen sich mit Hülfe der Gleichungen (38) ergiebt: 


2w 07 e ` ro 9 

— Yes — fe = 920,9) = l—20 +2g*—2g’+...)? 
= 2 w ARE 97 r N er [73 r9 a 
27)* SE Ve, — a = 0,9) ss OELäoLäatLäent. A 

TE 

2 w 


w 


-Ya = 3%,9) = Pati +gQ?+g°+....)] 


und somit nach (29): 


Ve ée De ek 
Zi nn: —— al SS 


LL ëSe £2g CL... 
29)* 2u 


Der unmittelbare Anblick der Gleichungen (29) und (29)* zeigt, dass x bezw. x’ in derselben 
Weise von q abhängen, wie x’ bezw. x von q’; umgekehrt ist also q dieselbe Funktion von x bezw. 
a, wie q’ von x’ bezw. x. 

Anstatt also, wie dies oben in den Gleichungen (30) und (31) geschehen ist, die Grössen A und 
q als Funktionen von x’ darzustellen, kann man auch A und o als Funktionen von x berechnen, 
indem man A definiert durch die Gleichung 

o 1 Lal 
50) A 2. —E; 
2 1-+Yx 


alsdann wird 
31)* EN = At ek 2 Ar 5 + 15 A3 E 150 an 13 十 1707 dii 17 L EE - 


Ist auf diese Weise der Wert von q’ gefunden, so ergiebt sich aus einer der Gleichungen (27)* 


+ 
e HI -~ ý Aw 
der Wert von SCH und hieraus nach Gleichung (35): 


r 
a 


AO e= S . Log t at 
q 
Die Reihenentwickelungen (31) und (31)* ergänzen sich hinsichtlich ihrer Konvergenz. Man 
wird die erste oder zweite dieser Reihen anwenden, je nachdem 
x > Vt =x, d. h, (da e + es + e= o ist) 


e So oder 


O > 
zZ A0 
SR 45 


E? 


e 7 @ o- igo: . ` INT A e S 
ist. Wenn daher 3 S 45° ist, so wird man der Reihe nach die Grössen 4 q und w mit Hülfe der 


e g É RAN a z 5 E 
Gleichungen (30), (31) und (27) berechnen; ist hingegen 3 > 45°, so wird man ebenso die Grössen 
t 
/ w e g - SÉ Do 
A, q’ und 一 aus den Gleichungen (30)*, (31)* und (27)*, und sodann den Wert von w mit Hülfe 
i 


der Gleichung (40) ermitteln. 
Damit ist auch die oben (pg. 11) aufgestellte Behauptung bewiesen: Für die meisten Werte 


von enügt bei Anwendung siebenstelliger Logarithmentafeln das 1% Glied der Reihenentwickelung 


DR 


e ; ? e a a ; a S 4 
(31) resp. (31)*; nur für verhältnismässig wenige Werte von — werden 2, niemals aber mehr Glieder 


jener Reihenentwickelung erforderlich sein. 

Die obigen Entwickelungen sind nun auch von Einfluss auf die Gleichung (26a), welche für eine 
gegebene Zeit t den Wert von up bestimmt. Die in dieser Gleichung auftretenden $-Funktionen können 
sowohl mittels der Gleichungen (17) durch eyklisch-trigonometrische Funktionen, als auch mittels der 
Gleichungen (38) durch hyperbolische Funktionen berechnet werden. Da die Reihenentwickelungen (17) 
und (38) nach Vielfachen von q bezw. q’ fortschreiten und daher um so schneller konvergieren, je 
kleiner q bezw. q’ ist, so wird man aus den oben angegebenen Gründen die Rechnung mit eyklisch- 


5 N = : k a : e ; 
trigonometrischen Funktionen vorziehen, wenn 3 < 45° ist, hingegen mit hyperbolischen Funktionen, 


Be gar z i : - 
wenn 5 > 45° ist. Es verdient noch hervorgehoben zu werden, dass die Rechnung mit hyperbolischen 
Funktionen leichter von statten geht, als die Rechnung mit ceyklischen Funktionen, da bei diesen auf 
den Wechsel der Vorzeichen in den verschiedenen Quadranten zu achten ist, was bei jenen Funktionen 
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i nicht der Fall ist. Dabei bliebe nur zu wünschen übrig, dass die Gronau’schen Tafeln, welche m. E, 
5? 
den Forti’schen („Tavole de’ numeri e delle funzioni circolari ed iperboliche“, Torino-Firenze-Milano 
》 ? 
1870) trotz der entgegengesetzten Meinung Günther’s bei weitem vorzuziehen sind, von D auf 7 Dezimalen 
vervollständigt würden. 


Numerische Rechnung. 

Ein Pendel von der Länge 1 = 1m besitze zu Anfang der Bewegung die Neigung oi < a < 7) 
gegen die abwärts gerichtete Vertikale und werde ohne Anfangsgeschwindigkeit eine Stunde lang in Be- 
wegung gelassen; die Beschleunigung der Schwerkraft sei g = 9,81. Wo befindet sich das Pendel genau 
nach Ablauf dieser Stunde? 

Hi 


1. Beispiel: Der Ausschlagswinkel « sei gleich 80°, also „ = 40°. Für diesen Wert wird 
x’ = cos 40° 一 0,7660444 
Val = 0,8752396, 
logA = 0,5219900 — 2; A = 0,03326519 
5log å = 0,6099500 — 8 


log2 = 0,8010300 
0,9109800 — 8: 23° — 0.00000008 


q = 0.03326527 
log q = 0,5219910 — 2 
J (o, OU = 1 + 2q 2 e 1,06653299 , 
mithin nach Gleichung (28a) 
log w 0,7562340 — 1; œ = 0,57047158 


T = 2» = 1,14094316. 

Die Dauer einer einfachen Schwingung beträgt somit 1,14094316 Sekunden. Um die Anzahl 
der Schwingungen zu berechnen, welche das Pendel innerhalb einer Stunde ausführt, hat man mit 
1,14094316 in 3600 zu dividieren und findet als Quotient 3155 mit dem Rest 0,82433. Das Pendel führt 
somit in einer Stunde 3155 einfache Schwingungen aus und fällt nachher noch 0,82433 Sekunden. 

Wo befindet sich nun das Pendel nach Verlauf dieser 0,32433 Sekunden? 

Da dieser Zeitraum kleiner ist als die Dauer einer halben einfachen Schwingung, so wird das 
Pendel aus der höchsten Lage noch so weit fallen, bis sein Neigungswinkel vu gegen die Vertikale 
den durch die Gleichung (26a) 


Lt 3 tz] = 5 ta 
| 3 6 
C08 - 4°cos | q” Cos u 
N, / Y(U) ` 1 20 ' | 2 w DES leg L 
sin = Var = Va, ZO: es e 
2 9 (u) i 2 ta Ve EE 
h 1 -+ 2 q cos Zu cos -< 
2 w 2 w 


bestimmten Wert erhält. Die Glieder hinter dem vertikalen Striche kommen bei Anwendung sieben- 
stelliger Logarithmentafeln nicht mehr in Betracht, in den meisten Fällen genügen sogar noch 
weniger Glieder. À 


` tor t. 180° ge 3 
Nun ist F = 5 STEE 186. == 5101073616 
w á w 
Lt SEET Mp eg 
cosg 一- sin 38° 49' 56 "4 = 0,62704357 
w 
2 tre 5 t Fr e 
Zu cos oo —- 2 q sin 12°20" 7”,2 三 一 001421311 
w 


| 


mithin, wenn man zur Abkürzung den Zähler des obigen Bruches mit Z, den Nenner mit N bezeichnet, 


q?sin 63030’ 10”,8 一 — 0,00099034 
vi sin 65° 19° 45”,6 一 — 0,00000223 , 


Z = 0,62605323 
N 0,98578466 . 
Da ferner x = sin = — sin 40° ist, so ergiebt sich 
3 log x = 0,9040337 一 1 
log 2 0,3010300 
i log q 0,6304978 — 1 
log Z = 0,7966113 — 1 
1 RECH 
log = 0,0062179 
= W ag H 
log sin 5 = 0,6383907 — 1 
y KOAR’ AAF AK 
5 25° 46’ 44,45 , 
und somit 
U = 51 0 33 28,9 Dr a 


2. Beispiel: Der Ausschlagswinkel 


In diesem Falle ist 


= 50:00 一 


also ebenso gross wie x’ im vorigen Beispiele. 


æ sei gleich 100°, also S = 


A 


cos 40°, 


Demnach wird auch 


A = 0,03326519 
o 0,03326527 ; log q = 0,5219910 — 2, 
und da zufolge der Gleichungen (27) und (27)*, nämlich 
Zu 
E Ve, P etc (1 2q 十 2q +...) 
zu x 
2 w e e 5 b 9 
gl Kg, Se J,” (0,9) (1+2q +2gqg* a 
7t 


’ 
w ` e 
- dieselbe Funktion 
i 


H 
w 


Zur Berechnung von 


w 


m 


0,57047158; 


von o ist, wie » von q: 


0,7562340 — 1. 


w 
log — 
Se 


w dient die Formel (40), nämlich 


r 


1 
mi Log S 


q 


Aus dem obigen Werte von log o ergiebt sich zunächst 


1 
log 
WE, 


1,4780090 


》 
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und hieraus mit Benutzung einer Tafel zur Verwandlung der gemeinen Logarithmen in natürliche 


1 = 
Log 一 = 3,4032415 
1 
1 ee 
log Log - ; = 0,5318927 
] 
log — 一 0,7562340 — 1 
i 
1 e EE 
log = 0,5028501 — ] 
-N 


logo = 0,7909768 — 1; œ = 0,61798343 . 
Die Dauer einer einfachen Schwingung beträgt also für den vorliegenden Fall 
dE 1,23596686 Sekunden. 
Daraus folgt, dass das Pendel innerhalb einer Stunde 2912 einfache Schwingungen ausführt und 
ausserdem noch weitere 0,8645 Sekunden in Bewegung ist. Der Ausschlagswinkel, den das Pendel 
nach Verlauf dieser Zeit besitzt, wird bestimmt durch die Gleichung (26a) 


j — A (u) t i 0.8645 
siny = yz -x Dee F ‚3645 
P ans TE 
und ist, da t grösser ist als der Wert der halben einfachen Schwingungsdauer, nämlich = w + 0,24652, 


negativ. Das Pendel wird demnach während der Zeit t den halben Schwingungsbogen zurücklegen 
und nach dem Durchgange durch die Vertikale noch 0,24552 Sekunden lang steigen. Bezeichnen wir 
den Winkel, den das Pendel während dieser Zeit beschreibt, mit 9, so ist nach Gleichung (23a) 


) F 4 
sin? = ee v=- , t = 024652. 
2 9 (v, q) E? 


Führen wir endlich mittels der Formeln (37) und (38) statt des Moduls q den oben berechneten 
Modul q’ ein, so erhalten wir 


n?2—=yY huet, ui Sin vir — q” Sin 3 var + oäBnmnbhvnr F... 
Sin 二 = Va: - 一 这 CNS 5 S ; z 
2 Ahv, q’) l Cos v'æ -+ q” Cos 3 v'm -+ q'® Cos 5 v'm 十 
7 Z s wi 
€ lt cia v 
N , w 
d rn 
also vm = ‚ Ist. 
2 w 


L 
. ? w . . 
Aus den obigen Werten für r und 一 ergiebt sich 
: i 


log Gezei 0,8317381 — 1 ; 
addiert man hierzu, um zu den Logarithmen der hyperbolischen Funktionen überzugehen, 
log M 0,6377843 — 1, 
wo M den Modul des Briggs’schen Logarithmensystems bedeutet, so erhält man 


log (ON vis = 0,4695224 — 1, 
also 


Mv’n 0,29479|7 一 0,29480 
3Mv'z = 0,88439 
5Mv 1,47398, 
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und hieraus mit Benutzung der Gronau’schen Tafeln: 
log Sin v'm 一 0,86460 —1 log Cos vn = 0,09320 


log Sin 3 v’n 0,57591 log Cos 5 v'æ = 0,59069 
log Sin 5 vin = 1,17246 log Cos 5 v'm = 1,17844. 


Da ferner 
log ou = 0,521991 — 2 


2. Jo oi 0,043982 一 3 
6.log d = 0,131946 一 9 
ist, so kommen in den obigen Reihenentwickelungen für Z und N wieder nur die beiden ersten Glieder 


in Betracht, und es wird 


Sin va 0,73215 Cos vr = 1,23937 
q*.Sin3 va 0,00417 q”. Cos 3 v'm = 0,00431, 
also 
Z = 0,72798 N — 1,24368. 
Setzt man diese Werte ein in die obige Formel sin Z UR und berücksichtigt, dass 
Z = SIN OOS, blogx — 0,94213 — 1 
ist, so erhält man 
Co 
7 


= 30° 49,1; also 9 = 61°39,2. 


Das Pendel führt also innerhalb einer Stunde 2912 einfache Schwingungen aus, fällt dann bis 
zur Vertikalen und beschreibt endlich noch einen Winkel von 61° 38°,2. 一 
Von praktischem Interesse sind insbesondere die Werte für die Schwingungsdauer des Pendels 
I gung 
bei verschiedenen Amplituden. Wie aus den beiden vorhergehenden Beispielen ersichtlich, braucht 
5 》 


。 a : ; a ` ; E 
man dieselben nur für Ausschlagswinkel bis zu 90° E 三 : 45°) direkt zu berechnen; die Werte der 


Schwingungsdauer für Ausschlagswinkel zwischen 90° und 180° erhält man aus jenen mittels des 
Überganges von den künstlichen zu den natürlichen Logarithmen. Dabei ergiebt sich, dass, während 


i a 2 8 e SE d w 
der Winkel 5 von 0° bis 90° wächst, œ die Werte von 0,50152 bis œ durchläuft. Für 5 = 45° 
u ei 
3 E ? ge ed = an N s: 
ist wœ = 0,59196, für 5 = HOT w 1,00679, für 3- 89° w = 1,73524; erst bei dem Uber- 
a a ` : ; . 
gange von „ = 89° zu e" 90° wächst mw von dem soeben angegebenen bis zu einem unendlich 


grossen Werte. — 


2. Fall ES > 21). 


Für diesen Fall gelten die Formeln (18) bis (20) und (13»), aus denen sich ergiebt 


4, 
sin s $ SI / LER 9, (u) u = I — to 
© —-@ 


J(u) ’ 2 w 
e 5 A 4 el 2 
4l) ai = = a Er e 
el Ga Vu" 


wo vo die Geschwindigkeit bezeichnet. mit welcher das Pendel die tiefste Lage passiert. 


3* 
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Nimmt man an, dass sich zur Zeit t = o 
dass der Winkel up gleichzeitig mit t wächst, 
das positive Vorzeichen zu geben; es wird demnach 


das Pendel in seiner tiefsten Lage befindet und ferner, 
so hat man tọ = o zu setzen und der Quadratwurzel 


D j , 
sin z - Í 1 Dé (u) i t 
= "x Hm)’ Ets 
Da si und x positiv und kleiner als 1 sind, so kann man auch in diesem Falle 
Nr; 
sin - 
S i "3 - Q 7t 
21») O < 
a 2 2 
x cosg 
setzen. Eine unmittelbare mechanische Bedeutung kommt hier dem Winkel w nicht zu: doch hängt 


von der Grösse desselben der Wert von v, ab, und umgekehrt entspricht jedem beliebig gegebenen 
Werte von v,, der grösser als 2 Vgl ist, ein eindeutig bestimmter Winkel æ. 
Aus (23b) folet mit Rücksicht auf die Gleichungen (17) und (27), dass 
? 5 d 


* Ui S 
iur Lo u 0, sn, = 0, also ur o wird; 
s Res A. (u) 9 (è) A (0 / . WY 
hingegen wird für t w u er ur = SCH Lë yx, sin“ l, also er. 
J (0) Da (È) 3; (0) 2 
Mithin ist die Dauer eines ganzen Umlaufs 
241) Jos Zo, 


Zur Berechnung von w kann zunächst wieder eine der Gleichungen (27) benutzt werden; die 


zweite derselben liefert, da 


V š 
€; Cy FRE ISt, 
lz zë 
w Fa (0, q), 
Vo 
und somit . 
na 2 | t 2 ] 元 
28b) 1 20 Seh Zo, 4) > +24 +2a +.. Sak 
0 Vo 


Den Wert von q erhält man wie im ersten Falle mittels der Gleichungen (30) und (31). 
Statt dessen kann man auch mittels der Gleichungen (30)* und (31)* erst die Grösse q’ be- 
rechnen, und sodann den aus (27)* resultierenden Wert von 


42) a 
1 


lz E e 
CN 


0 


Der erste Weg wird im allgemeinen vorzuziehen sein, wenn 


ə e a w ot 
x? $ = 45°, also Yo > 8gl 
ist; hingegen wird man den zweiten Weg einschlagen, wenn 
> 4 5°C e 080% also dgl ei < 
x č, D ER , also g h s 8 gl 
ist. — 
E EE ER i 全 全 21m 
) Für q = o geht diese Gleichung über in die bekannte Formel für die Centralbewegung: T = 
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Numerische Rechnung. 


Es werde wie im vorigen Falle 


L= Fm ud e =g 
angenommen; dann ist y3 gl 8,858894. Je nachdem also 


vo = 8,858894 
ist, wird man der Reihe nach die Grössen 


A A ® 
7 
w 
bezw. 4, og =, W 
i 


berechnen. 


Nehmen wir, um das Schellbach’sche Beispiel (1. c. $ 103) vorzuführen, an, es sei 
Vo dës 

d. h. der materielle Punkt besitze im tiefsten Punkte der Kreisbahn eine horizontale Anfangsgeschwindig- 

keit von 7 m. Wo wird sich alsdann das Pendel nach Verlauf von 1 Stunde befinden? 


Nach Formel (41) ist 


Cr 4.98 
Bi} S1. „— 1.3934 = 0,8948834150495 
vo” 49 í 
Ve 一 0,945982777354 
1 = IX 一 0.01387916|257 = q, 
á 14 yx 
also log # = log q = 0,1423633 — 2 


3,(0,9) =A 2 Lët iz 1,02775840 


und hieraus nach Formel (42) 
D 


log — = 0,6758339 — 1. 


1 
e z R `; ; SEN ] a I, 
Zur Berechnung von w dient wieder die Gleichung (40); dieselbe liefert, da en — 4,2773666 ist: 
log o 0,8098605 — 1; w 0,6454469 . 
Somit beträgt die Dauer eines ganzen Umlaufs 
T 一 1,2908938 Sekunden. 

Da nun 3600: T gleich 2788 mit dem Reste 0,98809 ist, so macht das Pendel innerhalb einer 
Stunde 2788 volle Umläufe und ist alsdann noch 0,98809 Sekunden in Bewegung. 

Während dieser 0,98809 Sekunden durchläuft das Pendel noch die halbe Kreisperipherie, gelangt 
also aus der tiefsten Lage in die höchste, wozu es 0,64545 Sekunden gebraucht, und bewegt sich 
dann wieder abwärts während eines Zeitraums von 

E = t— o - 0,34264 
Sekunden. Setzen wir daher in der Formel (23b) 


t d e w 
so wird, da 
J (v + è) - J (v) 
r +H) = lv) 
ist: 
— WY DC: 826 7 : 
sin * T EE , t = 0,34264. 


EA 
2 x Ze, ol 2 w 
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Auch hier empfiehlt sich wieder die Rechnung mit hyperbolischen Funktionen. Durch Einführung 
derselben mittels der Gleichungen (37) und (38) ergiebt sich 


了 / Ar H N r 
f S AG O IL 1 — 29’ Cos 2 vm + 2g4% Cos4vr|- ... 11 2 
sine — ere E (a L 1 — |2. 
2 x Jett, q) x 1 -+ 2q Cos? vm -+ 2q*Cos4vr |- ... E 
wo v GEN ist 
RE d, 
2 w 


' 
w 3 s ` ; z 
Setzt man für v und die oben angegebenen Werte ein, so ergiebt sich 


log (via = 0,0551241, also da 
logM — 0,6377843 一 1 ist, 
log (Mv’z) = 0,6929084 — 1, und hieraus: 
Mein = 0,49307 
2 Mvr = 0,98614 ; log Cos 2 v'm = 0,68972 
e A Mia = 1,97228; log Cos4 v’r = 167180 . 
Nun ist logq’ = 0,14236 — 2 und 4logq’ = 0,56945 — 8, mithin wird 
Zu Cos2v’n — 0,135866 
2 ou? Cos4v’n 0,000003481 ` 
also 
Z = 0,864137; log Z — 0,93658 — 1 
N = 1,135869; logN = 0,05533 


4log x 0,97588 — 1 
TL e 
log sin 5 — 0,90537 一 1 
un L 58988‘ 
2. 1280887. 


Von diesen beiden Werten ist der letztere zu nehmen, da, wie oben gezeigt wurde, das Pendel 
bei der letzten unvollständigen Schwingung noch mehr als einen Halbkreis zurücklegt. Es ist demnach 
W — 20256. 

Das Pendel hat somit, von der vertikal nach abwärts gerichteten Ruhelage ausgehend, während 
einer Stunde 2788 volle Umläufe vollendet und ausserdem noch 252° DÉI zurückgelegt. 

Das erhaltene Resultat weicht von dem Schellbach’schen (w — 252° 17’) nicht unerheblich ab. 
Man würde jedoch fehl gehen, wenn man diese Abweichung auf Konto der von mir angewandten 
hyperbolischen Funktionen setzen wollte; vielmehr ergiebt die Rechnung mit den gewöhnlichen sieben- 
stelligen logarithmisch-trigonometrischen Tafeln auf Grund der Formel (23%) bis auf die letzte Dezimal- 


stelle genau das oben erhaliene Resultat. — 


II. Abschnitt. 


Die Bewegung eines materiellen Punktes auf einer Kettenlinie. 


Die Gleichung der Kettenlinie in Bezug auf ein rechtwinkliges Koordinatensystem, dessen 
horizontale x-Axe vom Scheitelpunkt der Kurve den Abstand a besitzt, und dessen vertikal nach oben 
gerichtete z-Axe durch den Scheitelpunkt geht, besitzt die Form 


x 


a F Er e X 
2 =; Lo +e d a. Cos oder umgekehrt: 
E a 


I E | pa 1) 
Lra. Logl/ 3— II, 
TNA a? 


Das Prinzip der Erhaltung der lebendigen Kraft liefert die Gleichung 
} H g 


; e -A =): Sen 
SA Z ES d E SES DB, 


aus welcher sich, da nach Gleichung (1) 


1) 


lz x dx EN Ix 
ES is 9 and umgekehrt: 
dt a dt a” dt 
dE: 3 dz 
dt yz? er a” dt 
ist, ergiebt: 
ar N 3 d 
(; a? 1) E SE 
j ER 
d. h. wu 人 =—2g th, 
z°— a” \dt ; 
zdz dz 
3) dt = - EE EE $ 
y (z? — a?) (— 2 gz + h) YR.(z) 


wo R(z) das Polynom dritten Grades 
4) R(z) = (z?— a’) (— 2 gz + hb) 
bedeutet. 
Zur Bestimmung der Konstanten h nehmen wir an, dass im untersten Punkte der Kurve, also 
eE iee y S 
sei; dann ist 
—2ag +h, h =v? 
und somit 
R(z) = (z?— a?) . |v — 2g (z — all, 
Dieser Ausdruck verschwindet offenbar für 
上 
Bezeichnen wir diese letzte Wurzel zur Abkürzung mit b, setzen also 
9 


s Var 
5) a + — 


z EE 
2g $ 
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so wird 6) R(z) = 2g(z + a) (z — a) (b — z2) 
= — 2 gz’ + 2bgz* + 2a’gz — 2a'bg. 一 


z = Së P dz A . 

Es handelt sich nun zunächst wieder darum, den Differentialausdruck — —— auf die kanonische 
/及 (2) 

d ds d S das P F 

Form —— zu bringen, wo S das Polynom 4 (s 


yS 
im I. Abschnitt gezeigt wurde, am einfachsten mittels der Substitution 
ER , 1 - 3 R' (zo) 


ve +R” (z7); Z= % T- 
r s | a 0/3 0 i R 1 Hf, 》 
S Zo 24 Sr R” (zo) 


&,) (8 — e) (5 — e€) bedeutet. Dies geschieht, wie 


wo unter z, irgend eine Wurzel der Gleichung R(z) = o zu verstehen ist. Man wählt zweckmässig 
für z, einen Wurzelwert, den die Variable z im Laufe der Bewegung anzunehmen imstande ist, also 
entweder z, =a oder Zọ = b (nicht etwa , = —a, da die Variable z für kein reelles t diesen 
Wert anzunehmen vermag). Wählen wir z, = a, so wird 


~ Ria L R’ (a 
i) Dee ` EI + R” (a); 2 =a 二 - 1 oe H 
Sc 24 s—;R”(a) 
und es ist, wie sich durch Differentiation der Gleichung (6) ergiebt: 
| R’(a) = 4ag.(b—a) = 2 av” 


8 
| R”(a) = 4g.(b— 3a). 


: E des d Eë Esch . 
Die Substitutionsgleichung (7) führt den Differentialausdruck ES i über in rg? d. h. es ist 
y R(z 二 
dz ds 
gh 三 = +- 
Ra) ~S 


Hieraus folgt mit Benutzung der Gleichungen (T) und (6) 


de? L R’ (a r ' 
g= ( ) RG = I al -2g (z -+ a) (z — a) (b — z) 
Sab 


dz 
2g D'RE. Ve a) (b — 2) 


o 
z— a)? 


Bestimmt man diejenigen Werte von z, für welche dieser Ausdruck verschwindet, und berechnet 
mit Hülfe der Substitutionsgleichung (T) die zugehörigen Werte von s, so stellen diese letzteren die 
Grössen e, € und es dar. Nun verschwindet aber S für z = 一 a z= b und z = %; die ent- 


sprechenden Werte von 8, die wieder einstweilen mit e. e”, e” bezeichnet werden mögen, sind 


》 


eg 3.1 Beat 二 2 

= TER TER: ES 
In. vi a e 1 1 
d em EE ee 
di 1 rr 1 
一 zR (a) = e bg SH? 


ordnen und erhalten: 


1 1 
HR bg +- ag el 一 6 = ag 
1 1 1 g 
r — P d A7 U A 一 e ken eg 
10) ug g?g 538 gar (a) | el 一 88 2 (V a) 
i o S y 2 
一 Am D E E = H. 
ae bg e, es > (b — a) 4 


Nachdem somit die Wurzeln des Polynoms S bestimmt sind, führen wir noch eine Variable u ein 
mittels der Differentialgleichung 


ds 
11) == St D 
yS 
deren vollständiges Integral 
s = p (u — u) 


ist, wo u, eine willkürliche Konstante bezeichnet. Dann wird 


t R’(a) ag (b — a 
ole A} | s j e - a 十 2 (b -a) = 
pu—w) —;R (a) p (u — u) — e 


und somit nach Gleichung (3) 


Ä ao l mi 
3) dt = zdu = (a a), ) a. 


p (u — w) = es 


Die Gleichung (12) liefert z als Funktion von u, während mittels der Differentialgleichung (13) 
u als Funktion von t bestimmt ist. Nach Integration dieser Differentialgleiehung wird demnach auch z 
als mittelbare Funktion von t bekannt sein. 
Um das Differential z.du integrieren zu können, unterwerfen wir den obigen Ausdruck für z 
einer einfachen Transformation. Wir setzen zunächst die willkürliche Konstante 
ID ie 
Dann ist 


ag (b — a) 
Eé a 
人 一 9 
und dem Werte u = w entspricht der Minimalwert von z, nämlich z = a, 
dem Werte u = o hingegen der Maximalwert von z, nämlich z = b. Denn für u = o ist 
p (u — w) o (w) e; da nun nach (10) e — &, = ag ist, so wird in der That für u = o 


z = a  (b— a) b. 
Da auch im vorliegenden Probleme alle drei Grössen e, e und es reell sind, so können die 
Perioden 2 w und 2’ der Funktion ø (u) mittels der Gleichungen 


olw) = @ Plo -FH w) = ge, p(w’) = &,; 
f 


: w Së 2 d d 
stets so bestimmt werden, dass w und ~ positiv reell sind. Während also z von b bis a abnimmt, 
i 


d. h. während der materielle Punkt von seiner höchsten bis zur tiefsten Lage gelangt, wächst u 
von 0 bis w. 
Nach der Formelsammlung Art. (19) ist nun 
(et — €o) (el — gel 


on 1 w) — Du = PA) KS? 


15) 


(el 一 es) (es 一 es) 
pu 一 es 


| 


puwto)— e = 


Aus der ersten dieser Formeln folgt 


M | (es — ea) (e, — 6s) plu) — e, _ oul 一 es 
A e re pu) 一 e ae EE ein — e ’ 
mithin wird 
z = a + (b— a) We _ a-+(b—a)- (1 一 一 ) 
oul — e, p( 一 es 


und dieser Ausdruck geht infolge der zweiten der Gleichungen (15) über in 


z = a + (b— a): (\ Bir uo) Si E Ne, ea — p (u +w’) 


ea 一 65s 


2 RE 2 2 : 
a [es 一 bu 十 o)] = (a el 一 一 "pu 一 ww) 


g s £ 
de h: 
b 
16) z 5 一 u -+ ow’) 
2 g 
Folglich wird 
N 2 b,2 d duto 
dt EH ett ) | dp. = | = Pan, w | du. 
3 ge | |3 g du o w- w’) | 


Dieser Ausdruck ist ein vollständiges Differential. Setzen wir zur Bestimmung der Integrations- 


konstanten noch fest, dass 
für u 0, d: b für e bet 0 


ist, rechnen wir also die Zeit von einem derjenigen Momente an, in denen der materielle Punkt seine 
höchste Lage besitzt, so ist die Anfangsgeschwindigkeit desselben gleich Null, und es ergiebt sich 


i b 2 | d (u +w’) o'(w) | 
一 >u -+ d 
3 g Jet +o’) "s(o)|’ 
oder da nach Art. (18) der Formelsammlung 
d K vu? 
y nu j n'u + 9 Po Si 
colu -+ w) = gie Je oz (u) Glo je Se 
o(o) 
und somit 
€ d u ) 
(u w’) ja oe Ee 
` Sé . 3 
g (lu w) I w Zon 9 ( u ) 
vo ES 
e aR. 
o (u ! w) olw’) nu ] Fo 20) 
o (u -| w) o (w) w - 2 w A. sl 
2 w 
ist: 
Als 
z Di 2 2 7 了 2w 
(od Ce, Se St u + 
D oi go ` VA 1 
Fo 2w 
Hieraus ergiebt sich die Zeit, welche vergeht, während u von o bis w wächst, z also von b bis a 
abnimmt, indem man u = w setzt. Für u = w wird aber 


ER a : u al { 
Po Ces al: o und % CES Aa 5) 2 (0) ; 
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daher ist die Dauer einer einfachen Schwingung, d. h. eines Hin- oder Herganges: 
b sën 2b 4 
ey) L ja + 1) 一 w =n 
> gw 3 g 
Führen wir noch diesen Wert von T in den obigen Ausdruck für t ein, so ergiebt sich definitiv: 
Q’ u 
mg, ge Po LG 》 q) 
19) t= 9 a apesi = A 
L 
F w gw I ( q) 


Ze? 


Um endlich auch den Wert von z durch 9-Funktionen darzustellen, gehen wir aus von dem 


oben abgeleiteten Ausdrucke 


Peg p(u) 一 el 


plu) 一 es 
und machen Gebrauch von den Eormeln: 
plu) 一 e = SE | 
` | Formels. Art. (17), 
0, (u) 
(u) E Ben 


sowie von den bekannten Formeln (Art. 34 und 35), die den Übergang von den o-Funktionen zu den 
$-Funktionen vermitteln. Wir erhalten alsdann: 


E CT 
g o,” (u) ` J,” (0) Ja (>) e — 6o Ja G) 
ER ZS 人 2) i [77 do S ` s SE J Et ) SCH) ëlo, y ji Sen ZG g TITA Ces 
Wen "IO Q, 2 Bo. d Q 2 
g x 由 j i “0 5) 


oder endlich mit Rücksicht auf die Gleichungen (10): 


20) z = a + V2a (b — a). 


Die Formeln (18), (19) und (20) enthalten die vollständige Lösung des Problems, falls es nur 
darauf ankommt, irgend eine analytische Form für das Resultat aufzustellen. Doch eignen sich die- 
selben nicht unter allen Umständen zur numerischen Berechnung. Vielmehr wird man, um möglichst 
rasch konvergente Reihenentwickelungen zu erzielen, bei denen höchstens die beiden ersten Glieder 
in Betracht kommen, jene Formeln für den Fall, dass x > x’ (d. h. q > 0,043214 > q’) ist, einer 
Transformation unterwerfen müssen, bei welcher der Modul q durch o ersetzt wird, und demzufolge 
an Stelle der eyklischen Funktionen wieder die hyperbolischen Funktionen treten. Insbesondere gilt 
dies von der Gleichung (19), welehe in Beziehung auf u transcendent ist, und welche, ebenso wie die 
ihrer Gestalt nach einfachere Kepler’sche Gleichung bei der Bestimmung einer Planetenbahn, „propter 
arcus et sinus heterogeneiam“ (Kepler, Astronomia nova, seu Commentaria de stella Martis, Pragae 
1609, pg. 300) nur durch ein Näherungsverfahren aufgelöst werden kann¥*). Der in beiden Fällen 
(x <x und x > x’) einzuschlagende Gang der Rechnung wird am zweckmässigsten an folgenden 
Beispielen erläutert werden. 


9 (u) 
lu) 
entwickelung darzustellen, die nach den sinus oder cosinus der Vielfachen von t fortschreitet, haben bisher zu keinem brauchbaren 


*) Die Versuche, eine Gleichung von der Form t = «u + £- nach u aufzulösen, d, h. u durch eine Reihen- 


Resultate geführt. Denn mit der Aufstellung gewisser Integralformeln für die Koeffizienten einer solchen Reihenentwickelung 
ist für die practische Berechnung noch nicht viel gewonnen. (Of. Scheibner, Supplement zu oben genannter Abhandlung, 


Le pag. XXII). 
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Numerische Rechnung. 
Es werde angenommen: 


1) der Scheitelpunkt der Kettenlinie habe von der horizontalen x-Axe den Abstand 


a= #cm, und 
2) die Ordinate des materiellen Punktes zu Anfange der Bewegung sei entweder 
b = Lem (1. Beispiel) 
oder b 10 cm (2. Beispiel). 


Wo befindet sich, fragen wir nns wieder, der materielle Punkt genau nach Ablauf 
einer Stunde? 


1. Beispiel. (a 3, b De 


Aus den Gleichungen (10) ergiebt sich 


u Da — Ga b— a ] ʻa > 2 
x 一 = - z; Si l or z 
el — Ga b zra, bi 3 


Da somit x’ > x ist, so wird es sich, wie im vorigen Abschnitte gezeigt worden ist, hier empfehlen, 
der Reihe nach die Grössen 2. q und œw zu berechnen, und, da q < q’ ist, bei den folgenden trigono- 
metrischen Rechnungen die cyklischen Funktionen zu Grunde zu legen. 

Durch zweifaches Radizieren ergiebt sich: 


1 1— yx 


0,025319892556 , 
9 CS ? 
a i e yx 


und hieraus: 
7 1 


215 二 15 和 9 二. = 0,02531991337 log q = 0,4034622 — 2 
Sina = 1 + 2q + 2qt +... = 1,05064064875; log Ziel = 0,02145420 , 


mithin, da nach der zweiten der Gleichungen (27) im vorigen Abschnitte 


7T D 
21) © = -Jg (0), 
2 re = 
im vorliegenden Falle also 
TU D t a 
22) o = 9, (0) 3,°(0, 9) 
2g(b-- a) V än 
ist: 
log» = 0,8056632 — 1; w 0,639239 , 


wo g wie früher gleich 9,81 angenommen worden ist. 
Zur Berechnung von 7 dient die Formel 
oO í 
23) T m? i ] =I 330° ST bio 


g 90 1 Ba? + de] (Formels. Art. 35): 


da. nun 
27 q? 0,01730965 ; "1125 që 一 0,000000033 


Su? 0,00192329; 5q = 0,000000001 
und somit 


1 — 27 q? + 125q" E... = 0,98269038 
1 一 3q?+ 5q°Ẹ... = 0,99807671 


ist, so wird nach obiger Formel 
log 1 0,1027080 7 = 1,2668. 
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Um endlich die Schwingungsdauer T zu ermitteln, wenden wir die Formel (18) an und erhalten 


2b a 4 TT 
ee 0,4261592 ; nu = 0,5165340 , 
: g 
mithin 
T = 0,9426932 Sekunden. 一 


Aus diesem Werte von T ergiebt sich, dass der materielle Punkt im Verlaufe einer Stunde 
3818 einfache Schwingungen ausführt und ausserdem noch 0,79736 Sekunden lang in Bewegung ist. 

Um zu ermitteln, welche Lage der materielle Punkt nach Ablauf dieser Zeit besitzt, haben wir 
in der Formel (19) t = 0,797836 zu setzen und den zugehörigen Wert von u zu berechnen. Ist dies 
geschehen, so ergiebt sich mit Hülfe der Gleichung (20) sofort der gesuchte Wert von z, welcher die 
Lage des materiellen Punktes zur Zeit t in eindeutiger Weise bestimmt. 

Die Gleichung (19) besitzt, ausführlicher geschrieben, die Form 


Le 
. Um au A HUN 
e q sin — — SO) sin + 3q’sin— +... 
24) t Tu s E: oi w w 
A Sa" ge e ur A Aur, 
i 1 — 2q cos -+ 2 qf cos — 2q” cos — T.. 
wo wW Ki 
wofür auch geschrieben werden kann 2): 
mm 2 $ 4 
E Tu , Aa ( urz o 2 un q . dus q . bus 
25) t en ae re — sin 十 ; Sin 4 1 ; sin k sik 
w Tw 1 — q” w 1 一 4q w 1 一 9 w 一 9 w | 


Beide Gleichungen sind zur Berechnung von u in gleich hohem Grade brauchbar. Wir benutzen 
die erste derselben; aus ihr ergiebt sich 


DN u 
U 元 zum , 
K Si cl ai | 
qsın —za sın ites 
26) 2 tw Hr l w l KI 2 w i Har Z 
40 u -全 Ser E — E P ES? a 
e zÄ 7 m m 
T gT 1 Sons um TETE- um 一 I gT N 
(NE w 


wo Z und N leicht verständliche Abkürzungen bezeichnen. 

Beachten wir, dass q im Verhältnis zum ersten Gliede auf der rechten Seite dieser Gleichung 
eine kleine Grösse ist, so können wir, um den ersten Näherungswert u, von u zu erhalten, q nebst 
seinen höheren Potenzen vernachlässigen und erhalten 

2 w 
u, T Big 
mithin, da 
log œ = 0,8056632 — 1 
og T = 0,9743704 — 1 
ogt — 0,9016544 — 1 


logu, = 0,0339772; u, = 1,0813771. 
Diesen Wert von u, setzen wir auf der rechten Seite der Gleichung (26) für u ein und erhalten 
so als zweiten Näherungswert von u: 


9 = 
u7 2 u7 


e Deg 
q sın — 2g°sin E e 
2 w KE l w i er GË 
Zeck: Gef Ze, D SEET ee re 
4 g u7 (E zT N 
>- 1— 2qcos = + 2 qt cos —-= F & À 
w w 


*) Cf. Scheibner 1. c., Supplement pg. III. 
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Nun ergiebt sich aber aus den obigen Werten von q, u und w: 


UF _ 5,314518 一 304°,49944 一 304° 29’ 58” 
2 
q. sin Y 一 — 002086694; äu, sin == = — 0,0000007674 
w HI 
了 2 7 
gu, om “7 一 0,02868231; 2q‘. cos "2 — 0,0000003079 ， 
w a 
also e 
Z, = 一 0.02086617; N, = 0,97131738 
Sn Z, 
. 1 — — 0,05838242 
T N 0,058382 
und somit 
u, = 1,1397595. 


In gleicher Weise erhält man als dritten Näherungswert 


e Hei WR Nuet | 
` Q SIn - 20 Sın St S 
IT 2w Ha l w 1 a — e 8m Z 
eT. ST F ES "SC EEN 
el 1 — 2 q cos Ba j 2 qf cos sr S 3 
w 0) 
wo 
Tei 一 3200,989 一 320° 56’ 20,4 
w 
ist. Demnach wird 
. UW SE 。 4 Us7 
q . sin 2 一 — 0,01595528; 2 4 .sin 一 一 “一 一 0,0000008044 
w w 
197€ & 2 uz í 
2q.cos ff = 0,03932059; Zou, cos 一 0,0000001692 , 
w w 
also 
Z 一 — 0,01595448; N, = 0,96067958 
8 元 Z, 
—— . — 0,04513408 
gT N, 


und somit 
Hat ` 


— 1,1265112;; 


Us 
w 


317°,20701 = 317° 12’ 25”,24. 


Setzt man die Rechnung in dieser Weise fort, so gelangt man zu folgender Reihe von Näherungs- 


werten: 
Da u, 7e d 0 1 KQ 
u = 1,0813771 ; — 304° 29° 58 
w) 
aaorkak Ugam auf RE’ Ont 
Ha = 1,1897595 ; 一 320° 56’ 20 ‚4 
w 
Pe Us7 947 e Sne 
u, = 1,1265112; —- = 317° 12’ 25”,24 
w 
ETS 2180 107 19” 5 
u, 一 1,1299263 ; — 318° IO 127,54 
w 
DANLLE U-2t 29170 REI mo 
u, 一 1,1290663; —— = 317” 55 407,8 
w 


au 


ien "Mag State 
ge = 1,1292839; —- = 317° 59’ 21” 
w 
(E ; 3 Bern 
gn, = 1,1292291; —— = 317° 58’ 25”,7 
w 
EN Us 和 ae d 
us = 1,1292429; == = 317° 58’ 40’ 
w 
: Wat EE 
u = 1,1292393; = 317° 58’ 36” 
w 
(TE: EE R S 
up — 1,1292403; 一 317° 58^ 37” 
w) 
or Uj 7 2170 ar 997 
u, = 1,1292401; 一 一 317° 58’ 37 
w 
i ' U. s = TATE 
Da dieser Näherungswert von u,, bezw. —!— mit dem folgenden übereinstimmt, so stellt derselbe 
w 
denjenigen Wert von u dar, welcher der Gleichung (24) Genüge leistet, und es ist somit 
uyr =0 OI DrA/ 
u = 1,1292401; — = 317° 58’ 37”. 
wW 
Endlich liefert die Formel (20): 
' e 2 
d ur ~ Aug 
| Zut coss 2 q4 cos SE 
2 w 2 w 
z = a 4+ Yalb—a)- Ju 
i f ur u —_ 
| 1 — 2a cos "7 H 2 qt cos 一 一 pidi | 
w w 
für den zuletzt angegebenen Wert von u ist aber 
2 ua Gef ` 3 um f , 
2 qa c08 = — 0,7447557; Zut cos 一 — 0,000231927 
Zo 2 w 
S ur EEN GU d 
2 q cos 0,0376191; Zu cos — 0,000000085 ， 
w) w 
mithin 
1 e H 7 
17t Ə UT = 
244 coss = 十 2q4 cos — 0,7449876 
uye ; 2 ur 
1 — 2 q cos 2 qf cos er 0,9623810 
w) w 
und demgemäss, da a = 0,5 und 2a (b a) = 0,5 iet: 


z = 0,92373. — 
Die Bewegung geht also in der Weise vor sich, dass der materielle Punkt während einer Stunde 


rg 
3818 einfache Schwingungen ausführt, sodann, da t >, ist, bis zum tiefsten Punkte der Bahn fällt 
und von da aus endlich noch wieder steigt, bis er den Abstand z = 0,92373 cm von der Abseissen- 
axe erreicht. 

2. Beispiel. (a = A b = 10.) 


a) Ableitung der zur numerischen Berechnuns geeigneten Formeln. 
Für diesen Fall wird zufolge der Gleichungen (10) 
19 S 2 


— =; WI LN 
e Bag e 21 ， 


ee 一 6 b— a 


x z 
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also umgekehrt wie für das vorige Beispiel: x > x. Aus den im vorigen Abschnitte dargelegten 
Gründen werden wir daher hier der Reihe nach zunächst die Grössen A, q w und alsdann w berechnen. 
Zur Berechnung von 4, q’ und w dienen die Formeln (30)*, (31)* und (40) des vorigen Ab- 


£ 
N S $ N) S ` GC 7 
schnittes, während sich zur Berechnung von 一 aus der zweiten der Gleichungen (27)* des vorigen 
q 3 g 
Abschnittes die Formel ergiebt 
r 
四 HI H e Hi 9 
27) — = ——— J (0, q) = 7 一 一 一 一、 9°(0,4Q). 
G K ? < s KR 
1 2 Ve, 一 es 2g(b 十 al) 


Die Ermittelung der Schwingungsdauer T erfordert ausserdem noch die Kenntnis der Grösse 3 
Zu dem Zwecke berechnen wir zunächst die Grösse 7’, welche von q’ bezw. w in derselben Weise 
abhängt, wie 7 von q bezw. œ, und welche daher bestimmt ist durch die Gleichung 


ai 1—21? +1254 F... 


T = ech Fre ` 
l 人 
oder 
28) vi EE: WE KT E n.. 
á ) 1 e S 3 = 76 H 
1 Ken E ES 


sodann ergiebt sich n mittels der bekannten Relation: 
zei 


29) no’ — no = Sep 


Ist auf diese Weise T gefunden, so hat man noch mit Hülfe der Gleichungen (19) und (20) den 
Wert von u und endlich denjenigen von z zu bestimmen. Führt man in die eben genannten Gleichungen 
statt des Moduls q den Modul q’ ein, so treten daselbst statt der eyklischen Funktionen wieder die 
hyperbolischen Funktionen auf. Aus den Gleichungen (37) und (38) des vorigen Abschnittes ergiebt 

i 


a INPE a d w { 
sich nämlich, wenn für den Augenblick — » v = v’ gesetzt wird: 
z w 
wni 。 RER 
beer éi Ae VE e 
fe = ne w . vg o C el £ 
dng, OU = è Ze, J (vi, q) 
also ; 
9 
I (v, q) 2 wri | wmi 2q% Sin vr Goar Sin 3 v'm +... 
小 (v q) 17; w W g 1 9 2 
dk, 2q’4 Gosviar 十 2q 生 Cos3vr +... 
d bh, 
u . umi OUT SE LE AW 
A (: 可 . Sina; + 3q” Sin — + 5q! Sn —— +... 
2 w umi , wi 2 w 2 w 2 w 
A S ) O e o uC c Dn 
( ‚O8 = 07 VO8 = f H 08 é 
0 2o 20 a SÉ ; EIN 
mithin geht Gleichung (19) über in: 
g; Umi a va 3 uri | 
D K Ə 5 D I ir Ji Dës 
Gr = zei N SLA 2 | EEN 
: 2 geil w ge uri d 3 um 
= = Cos H q” Cos 二 一 - 
2 w 2 w 
oder 3 
w Ui E e 
Sin CES 87 E oai 7 
w 2 / 
30) t = «u+ $ 一 z cu 十 6 去， 
uri 3 uri N 


dk Ke 
& = S 
2 w gww’ 
31) S 
m 
= BE 
gw 
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Die Gleichung (30) entspricht der obigen Gleichung (24) für den vorigen Fall, wo x < x’ war. 
Ihr zur Seite wollen wir nun noch eine andere Gleichung stellen, welche zu (30) in derselben Be- 
ziehung steht, wie Gleichung (25) zu (24), und welche für den vorliegenden Fall zur 3erechnung 
von u geeigneter ist, als die Gleichung (30). 

Aus der Theorie der elliptischen Funktionen ist bekannt*), dass sich der logarithmische Differential- 
quotient einer $-Funktion in doppelter Weise darstellen lässt, nämlich entweder ale Quotient zweier 
Fourier’'schen Reihen oder als einfache Fourier’sche Reihe; so ist insbesondere 


1 (vg) ` sin væ + 3q? sin äer + .... 
m A.C, ol co8V75 十 o eoeäer + 
q” i 
= —tgvyn SÉ ; sin 2 vr — Bin. o aa a eoo N 
S Lo Lag — 
also 
ah e GEN: dä € e i 2 4 
sin væ + 3g” sin äert, Dh s € 5 
L : - tg vor — 4 d | a Sin Zen — 1 4 Sin 4 vze 4 5) i 
cosvæ-}- q*cos3 vm +... 1—q l— g z 
Ersetzen wir in dieser Gleichung q durch a’ und v durch vi, so erhalten wir, da 
54 | D 
cos (vi) Cosv | e SE 
Ee G also tg (vi) = i Tg vy 
sin (vl) = i Sin v | 
ist: 
Sin væ +3g”Sin3vr +... HS EE e i 
-一 一 一 --， CH — = Teva —4- l z Sin 2 ve 一 a 77 Sin var + ..); 
Cosva+ qg’Cos3vr +... I — q? UG] = 


zufolge dieser Relation nimmt die Gleichung (30) die Form an: 


: P u S eh 
f uci 4q Umi E erhalt 
22) t = au + o (re > — > Sin 一 一 +- 1 = Sin —— F.J, 
E e I 2 1 TRN 4 
á UI q w ( w 


welche sich für den vorliegenden Fall zur numerischen Berechnung von u als besonders geeignet 
erweisen wird, 

Was endlich den Wert von z anbelangt, so ergiebt sich aus Gleichung (20) mit Rücksicht auf 
die Formeln (37) und (38) des vorigen Abschnitts: 


i uzri 2 uzi 5 
1 1— So Cos 4 Sen Ge ZA. 
[477 w 

29 Er A CR d 

33) z= a+ 2a b— a 5 l ge 

l Ayo’ ui SE 5 uzi 
tyq Cos ; + q” Cos + q’? Cos zk, 
2 w 2 w 2 w 


b) Numerische Berechnung, 
Wie bereits oben angegeben, ist für unser Beispiel 


19 = 2 
und wl = = 


92 PIE 


9 


x 


*) Cf. Scheibner, l. c. pg. IIT. 
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also Vx = 1,97528956275116 
É 1—}Y SE 
KES 1 ee 0,006254880844 
2 1-yx 
Pa a 0,000000000019, also nach Gleichung (31)* des vorigen Abschnitts 
a’ = 0,006254889863 ; log vu = 0,7962197 — 3. 
Hieraus folgt 
2q' = 0,012509779726 
2 q4 一 0,000000003061 
0,9) =1+27+2g° +... 1,012509782737 ; log 9, (0, q) 0,00539923 . 
Nach Gleichung (27) ist aber 
= F P k k I (0, q) , 
1 y2g (b-a) 
mithin wird, da 2g (b -+ a) = 21 g ist: 
log — = 0,3510042 —1; = = 0,2243904 . 
e 1 
Ferner ist nach Gleichung (40) des vorigen Abschnittes 
d ] 
wW = = ` Log ns 
gu 1 
nun ist 
1 


log o = 0,7962193 — 3; log = 2,2037807 , 
: d 
also 
l er 
Log : = 50743926; log Log 一 = 0,7053841 
4 b 1 
und somit i 
log œ = 0,5592384—1; œ = 0,8624419 . 
Zur Berechnung von ai und » dienen die Formeln (28) und (29). Aus Gleichung (28) ergiebt sich, da 
27q” 0,0010563385 
0,00011737094 
ist, während die Glieder mit q’® auf die zehn (bezw. 11) ersten Dezimalstellen ohne Einfluss sind: 
1 — 2794? + ... = 0,9989436615 


D 7/9 
ə] m 


1 一 3q? +... = 0,99988262906 , 
folglich 
1—27 q” ang 
d 一 0,999080922 , 
1 一 3q! 
mithin 


log (ni) = 0,5637063 ; 
hieraus folgt mit Hülfe von (29) 
log 三 0,0356215; ı 1,085479 . 


Setzen wir endlich die soeben gefundenen Werte von œ und n in die Gleichung (18) ein, so 


g = 


ergiebt sich für die Schwingungsdauer T der Wert 

T 一 2,85888 Sekunden. 
600 
T 


Punkt im Verlaufe einer Stunde 1259 einfache Schwingungen aus und befindet sich nachher noch 


3 S S bé = ; ; : ; 
Da nun 1259 mit dem Reste 0,67 ist, so führt im vorliegenden Falle der materielle 


0,67 Sekunden lang in Bewegung. 


Unsere nächste Aufgabe besteht darin, den dem Werte 
t = 061 


entsprechenden Wert von u mittels der Gleichung 


` E Ee 49, E ET EE 
32) e = u P a Tg ES e Kies T3 D 7 1 [71 Sir D I LK 
2 w 1 一 4 wW 1 一 4q w 


zu bestimmen; die konstanten Koeffizienten auf der rechten Seite dieser Gleichung besitzen die Werte 


dk 7 4 e ` 一 P 
ge Gs es 0,0062546 ; log o 一 0,7961995 — 3 
2 IO Boom e 
Ce mi. 1,4271734; log d 一 0,1544767 
on 
34) 6° d r 
d — 0,0000391251777; log-— 1; = 0,5924564— 5 
l— q“ li 
D WC 
7 0,00000000153066 ; log e 0,1848787 一 9 
1—q i a 


Wegen der Kleinheit der beiden letzten Koeffizienten werden wir dieselben zunächst vernachlässigen 
können und einen schon ziemlich genauen Näherungswert von u erhalten, wenn wir 
DE uw 
35) t=oeu+Pß.Tg 
setzen. Aber auch diese Gleichung ist in Beziehung auf u transcendent. Beachten wir jedoch, dass 
œ im Verhältnis zu ĝ eine kleine Grösse ist, so werden wir den ersten Näherungswert von u, erhalten, 
indem wir œ = o setzen; alsdann wird 
er SS UECHT t 
36) De 过 Kent 
2 w Héi 
also 
, Hut) 


log Tg ni 0,6715981 — 1 = 0,67160 — 1. 
2 w 


£ 


Dazu gehört nach den Gronau’schen Tafeln das Argument 


ui 
M . 2; = 0,22122, 
2 w 
mithin ist 
log u, = 0,8619244 —2; u, 0,0727653 . 
Aus Gleichung (35) folgt nun 
uyi t— ou 
37) Tg: : - : 
) oc w d 


wir erhalten demnach zur Bestimmung des zweiten Näherungswertes u, von u die Gleichung 
üri t— eu, 
38 Tee 一 e 
) I d 


dieselbe liefert, da 


t— au, 一 0,66954488 ; log (t— eu) = 0,8257797 — 1 


ist: 
log Tg 5, 一 0,6713080—1 = 0,67130—1, 
á w 
M-27 — 0,22104, 
á w 


log u, = 0,8615709— 2: u, = 0,072706. 


36 
Der dritte Näherungswert, den uns die Gleichung (37) liefern würde, stimmt mit u, überein. 
Wir werden daher zur Ermittelung des dritten Näherungswertes nunmehr die exakte Gleichung (32) 
benützen, aus der sich ergiebt: 


e Wi ec Ki O pe 
i url t— eu : ( . um q ~. UE) 
39) Tg5—; = 一 十 4 ( | LE Sin te A : 
2 w ß 1— q” w 1 一 d 全 Ce 
folglich haben wir zur Ermittelung des dritten Näherungswertes u, von u die Gleichung: 
(E t 一 au; 人 E 
m. Ms _ 2 | a, Ha 1 x GAN 
40) 1857 = zu + z Din -7 Fre 74 NIN n Zen . 
Í q w D w 


Aus dem oben gefundenen Werte von u, ergiebt sich aber 


t — au, g 
z = 0,4691408 
P 
M-27 一 0,4420794 — 0,44208; log Sin ZT — 0,080282 
Oo) d w 
2 ui en dë Zu GE 
M. -一 一 一 08841588 = 0,88416; log Sin 一 一 一 = 0,57567, 
0) OD) 
also mit Benützung der in (34) angegebenen Werte: 
4q”? oTt SCH 
A i= Sin — 0000188277 
zl on w 
Au? 2 uzei re 
I .Sin — 0,000000083 , 
1— q ` [0] g 
und demgemäss 
Tg 27 — 0,4693291; log Tg = 0,67148 — 1 
"Zo er 
u,7r1 
MN, 二 一 0,22115 
2 w 


log u, — 0,8617869 — 2; u, = 0,0727423 . 


Berechnet man in derselben Weise den vierten Näherungswert u, von u, indem man auf der 
rechten Seite der Gleichung (39) u durch den eben gefundenen Wert u, ersetzt, so erhält man 


Dutt I Us 和 fl 
Tg 2— 0,4693291 — Te = 一 
` SC eme 2 w 
es ist also bereits u, = u,, und somit ist derjenige Wert von u, welcher der Gleichung (32) Genüge 


leistet und dem Werte t = 0,67 entspricht: 
u = 0,0721423; logu =. 0,8617869—2 . — 
Endlich ist dieser Wert von u noch in Gleichung (33) einzusetzen und der zugehörige Wert 
von z zu berechnen. Für den soeben angegebenen Wert von u ist aber: 


VI = 0,44230; log Cos = = 0,19451; 2q Cos F = 0,0195777 
w w w 
M-== — 0,22115; log Cos == — 0,054029; Cos ZU = 1,13248 
2 w P > 2 w i 2 w 
3 uri Da 3 uci In Sum , 
M- EE 0,66345; log Cos we 0,38241 ; q? Cos -5—7 = 0,0000944 , 
IO w á w 


während die Entwickelungsglieder mit höheren Potenzen von q’ nicht mehr iin Betracht kommen. 


Demnach wird, da a = 0,5 und 2a.(b—a) = 95 
ist: 
2 = 09-1 7,50217 d.h. 
z =. SOME: 
Im vorliegenden Falle beschreibt also der materielle Punkt innerhalb einer Stunde 1259 einfache 
Schwingungen, gelangt nach Zurücklegung der letzten vollen Schwingung zu dem der Anfangslage d 


entgegengesetzten höchsten Punkte der Bahn und fällt von hier aus noch bis zu demjenigen auf der- 
selben Seite der z-Axe liegenden Punkte der Bahnkurve, welcher von der Abseissenaxe den Abstand 
7,8011 cm besitzt. 


In ganz analoger Weise, wie die Bewegung eines materiellen Punktes auf einer Kettenlinie, lässt 
sich auch die Bewegung eines materiellen Punktes auf einer Parabel bestimmen, deren Axe 
vertikal nach oben gerichtet ist. Nehmen wir die Axe der Parabel zur z-Axe, die durch den 
Scheitelpunkt der Parabel gelegte Horizontale zur x-Axe und hezeichnen den Parameter mit 2a, so 
lautet die Gleichung der Parabel: 

x“ 一 437 
und die Differentialgleichung zur Bestimmung von z als Funktion von t: 
dt = (z +- a). GE S 
VRG) 
wo 


e? . 
R(z) z (z +- a) (— 2gz + h) = 2gz(z + a) (4 al = 2gz (2 + a) (b—z) 
r ER > 


ist; vọ und b haben dieselbe Bedeutung wie oben: v, stellt die Geschwindigkeit dar, welche der 
9 

den Maximalwert 
2 D 


已 


materielle Punkt in seiner tiefsten Lage, hier also für z 一 o, besitzt, und b = ;° 


von z, der während der Bewegung erreicht wird. 

Da sich die Rechnung in ihrem weiterem Verlaufe genau so gestaltet, wie für das zuletzt be- 
handelte Problem, so wird es genügen, das Resultat derselben anzugeben. Die zur Berechnung der 
Perioden erforderlichen Grössen el pe, e€, besitzen nebst ihren Differenzen die Werte 


ag bg ag 
el = 一 ~ el 一 es = x 
3 6 F 2 
ag be a 十 hb 
SE — Ba 一 or 
2 6 | 6 Di €g 2 5 
ag bg bg 
e = — -F — | &— e = —; 
Di 6 3 | d bi 2 H 


endlich treten an Stelle der Gleichungen (18), (19) und (20) die folgenden: 


Auf Grund dieser Formeln vollzieht sich die numerische Rechnung genau nach den oben am- 


gegebenen Vorschriften. — 

Wesentlich anders gestaltet sich jedoch die Rechnung für die Bewegung eines materiellen 
Punktes auf einer Parabel mit horizontaler Axe. Bei diesem Probleme besitzt das Polynom 
dritten Grades R(z), dessen Wurzeln in den bisher behandelten Fällen sämtlich reell waren, nur eine 
reelle und zwei konjugiert imaginäre Wurzeln: z = 2ai und z = — 2ai; infolge dessen wird auch 


von den Grössen ej, gx und e, nur die eine (e,) reell, während die beiden andern (e, und e) kon- 
jugiert komplexe Werte erhalten. Die bisher angewandten Formeln bedürfen daher, wenn die Resultate 
in möglichst reeller Form erscheinen sollen, für diesen Fall einer durchgreifenden Transformation, 
deren Grundzüge im Artikel (46) der „Formeln und Lehrsätze zum Gebrauche der elliptischen Funktionen“ 
angegeben sind. Es war ursprünglich meine Absicht, auch noch diesen Fall der Bewegung in ebenso 
erschöpfender Weise an dieser Stelle zu behandeln, wie die vorhergehenden. Doch muss ich aus 
Rücksicht auf den dieser Programmabhandlung zugemessenen Raum davon Abstand nehmen, und ich 
behalte mir daher, da eine blosse Mitteilung der resultierenden und zur numerischen Berechnung ge- 
eigneten Formeln ziemlich zwecklos wäre, vor, bei einer späteren Gelegenheit auf dieses ebenso 
interessante wie lehrreiche Problem zurückzukommen. 


